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Abstract 

We study analytically the medium-induced energy loss of a high energy parton 
passing through a finite size QCD plasma, which is expanding longitudinally 
according to Bjorken's model. We extend the BDMPS formalism already 
applied to static media to the case of a quark which hits successive layers of 
matter of decreasing temperature, and we show that the resulting radiative 
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energy loss can be as large as 6 times the corresponding one in a static plasma 
at the reference temperature T = T(L), which is reached after the quark 
propagates a distance L. 

I. INTRODUCTION 

Recent work fl|-[8[ on medium-stimulated gluon radiation from fast partons traversing 
(hot and cold) QCD matter starts from the assumption, that the properties of the medium 
and its interactions with the energetic quark or gluon projectile do not change with time, 
i.e. the basic parameters /i, which is the typical transverse momentum given to the parton 
by a single scattering in the medium, and the parton's mean free path A are kept constant 
in time. It also means in particular that the temperature T remains constant during the 
time the parton is passing through the QCD plasma. 

In this paper we study analytically the propagation of a quark, of high energy E, travers- 
ing an expanding hot QCD medium, i.e. we investigate jet broadening, induced gluon radi- 
ation and the resulting radiative energy loss of the quark. Thereby we extend the analysis 
of Q to the case of time-dependent parameters [i and A. We follow BDMPS HH, and we 
take into account our recent work @], in which we also show the equivalence of our approach 
with B. Zakharov's |7|J| formulation of the Landau-Pomeranchuk-Migdal effect || for QCD. 

For simplicity we consider a high energy quark entering and passing through a hot QCD 
medium. We may imagine the medium to be a quark-gluon plasma produced in a relativistic 
central A — A collision, which occurs at r = 0. At time tq the quark enters the homogeneous 
plasma at high temperature T , which expands longitudinally with respect to the collision 
axis. We may consider tq to be the thermalization time. For most of our results the limit 
t —>■ can be taken with impunity. We shall also state results for the more realistic situation 
where the quark is produced at tq = in the (not yet thermalized) medium. The quark, for 
simplicity, is assumed to propagate in the transverse direction with vanishing longitudinal 
momentum, such that its energy is equal to its transverse momentum. On its way through 
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the plasma the quark hits layers of matter which are cooled down due to the longitudinal 
expansion. We also assume that the plasma lives long enough so that the quark is able to 
propagate a given distance L within the quark-gluon phase of matter. 

The properties of the expanding plasma are described by the hydrodynamical model 
proposed by Bjorken (TtJ. The parameters \i and A depend on temperature, and therefore 
on time. The main relation is the scaling law 

TV = const, (1) 

where r is the proper time of the expanding medium; at rapidity y — it coincides with 
the distance (time) the quark has propagated through the plasma. The power a, which we 
approximate by a constant, may take values between a = and a = 1 for an ideal fluid. 

Let us state our main results for an expanding medium with a < 1. As for the static 
medium the transverse momentum broadening of the jet follows the random walk behaviour, 
namely the characteristic width p\ w is proportional to the path length L. The radiative 
energy loss per unit distance, —dE/dz, can be as large as 6 (2) times the corresponding loss 
in a static plasma at temperature T = T(L). The number 6 (2) corresponds to the situation 
where the quark enters the expanding plasma from outside (is produced inside the plasma). 

This paper is organized as follows. In Section II we treat jet broadening due to multiple 
scattering in case of an expanding plasma (with a < 1), and we estimate the characteristic 
width p\ w . Section III deals with the induced gluon radiation. In Section IV we derive 
the energy loss of a quark and relate it to p\y/. Following Bjorken |10| we review the 



main characteristics of an expanding plasma in Appendix A. The Green function of the 
Schrodinger-like equation with the time dependent "potential" is studied in Appendix B. 
Integrals which are necessary in calculating the energy loss are presented in Appendix C. 

II. JET BROADENING IN AN EXPANDING MEDIUM 



In this section we consider a high energy parton propagating through an expanding QCD 
medium. By multiple scattering a transverse momentum is given to the parton. In || we 
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have summarized the derivation of the resulting transverse momentum broadening for the 
case of a static uniform medium. In the following we generalize this derivation taking into 
account the space-time development of the medium. As described in Appendix A we assume 
longitudinal expansion. 

Because of the evolution of the medium the parton propagating in the transverse direc- 
tion, z, is affected by the position-dependent density of the plasma p(z) and the parton cross 
section der/d 2 q±(q±,z). 

Based on the probabilistic interpretationF] the master equation for the probability f(q±, z) 
for a quark to have transverse momentum q± (orthogonal to its direction of motion) at 
position z is 

df(q 2 ±,z) r 2 da f _ , 2w 



dz j f{q±, z )p( z )j2gi(q± - 9±> z ) d <?j_' 

+ / f(q'l,z)p(z)^(ql-q ± ,z)d 2 ql. (2) 

The first term (loss term) accounts for partons which are scattered out of the direction 
q±,q± — > qj_i and the second one (gain term) counts those partons which are scattered 
into the direction q± from all other directions q_[, q_[ — > q±. The result given in is 
reproduced with l/ada/d 2 q± and the mean free path A = 1/ per independent of z , where 
a = J d 2 q± dcr/d 2 q±. With a ^-dependent mean free path 

X(z) = [p(z)a(z)]-\ (3) 

(0) can be written as 

m dJ %^ = -n&z) + f±^tfLz)M&-*L) 2 ,*)*& (4) 

which can be diagonalized by defining 



x The main difference from the static case is the expression for the absorption of the parton along 



its path between points zq and z: exp[— (z — zq)/X] for the static and exp —J* dz' p(z')a(z') 
the expanding plasma, respectively. 



for 



and 



f(b 2 ,z) = Jd 2 q ± e- rb ^f(ql,z), (5) 



V{b\z) = I d 2 q ± e-^J^-(q ± ,z). (6) 



The resulting equation becomes 



m? l ^ = -\#v(b\z)f(b 2 ,z), (7) 



where 



v{b\z) = Ul-V{b\z)l (8) 
b 2 

and V(0, z) = 1. As discussed in in QCD v(b 2 , z) has no finite limit for b 2 — > 0, never- 
theless, (0) may be solved in a logarithmic approximation 

v(b 2 ,z) ~ /i 2 (z)5, (9) 

independent of 6. As in |2|,[| we introduce the scale /i 2 , with /^(z) representing a typical 
momentum transfer to the parton in a parton-medium collision, evaluated at position z. An 
explicit model for the scattering cross section is given by the screened "potential" 

= ~4^= r^L 2,2 - ( 10 ) 

a d z q± 7i (q^ + fx 2 ) 2 

For b 2 ~ we get, using (|) 

with the (transport) coefficient || defined by 

The solution of ([TT|) is 

/(^,z) = /o(^zb)exp[-j /Vg(z')}, (13) 



from which the characteristic width of the distribution f(q±,z) is deduced 

Plw(z) = <£(*)> = f Z dz'q(z')- (14) 
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For a hot (massless) medium the z dependence of q(z) may be determined from the tempera- 
ture dependence of the expanding fluid, T = T(z). The leading term of the high-temperature 
expansion for q(z) in fll2|) is determined by the T-dependence of the density p(z) of the 
medium, 

q(z) = q(zo)(T/T ) 3 . (15) 
This implies that the medium undergoes cooling from T to T when the parton propagates 



from z to z. Using Bjorken's model [1(J summarized in Appendix A, we may write (see 
eq.M) 



. a 

^0 



q(z) = q(z )^) . (16) 

Let us consider the interesting case of an interacting and expanding plasma, i.e. the case 
a < 1. Inserting (|T6| ) into ( Tjj ) the integration (with a = const) gives for z = L in the limit 



z — > the random walk behaviour 



In general the relationship is 



pV(L) = g(L)L UJ , (18) 
1 — a 

which shows the delicacy of taking the limits zq — > 0, a — > 1. 

In the high temperature phase of QCD matter, we note that a = 1/(1 + Ax/3) and, since 
Ai = 0(af) (cf. eq.QATD in Appendix A), 

1 - a ~ Ai/3 = 0{a 2 s ). (19) 
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III. GLUON RADIATION SPECTRUM IN AN EXPANDING MEDIUM 



Here we generalize the derivation of the soft gluon emission spectrum M to the case of 
an expanding hot medium. As described in the Introduction we assume that the fast quark 
is produced by a hard collision outside the medium. Let us first start with the key equations 
- valid in the static case - of sect. 4 in 0, which are re-examined in [|]]. 

Because of the Landau- Pomeranchuk-Migdal phenomenon ||, the induced spectrum is 
determined by an interference, essentially by the gluon emission amplitude at ti, f(b, t 2 — t 1 ), 
evolved in time to t 2 > h, and the complex conjugate Born amplitude f *(b) for emission at 
t 2 . We keep all the variables unsealed, as we did in the previous section. The Born 6-space 
amplitude for gluon emission is given by 

— * 

fo(b) = -4^(1 - V(b 2 ))^ ~ -m/M (20) 

where we work in the logarithmic approximation (cf. eq.@). The two-dimensional vector 
structure of fo and / takes into account the two polarizations of the emitted gluon. 

The induced gluon radiation spectrum (per unit length), in the limit of soft gluon energy 
u and in the large N c limit (cf. eq.(4.24) in ||), is given by 

udl a s N c l j [L dt 2 r** dh r d 2 b -,-> r*f^ K= °\ (o^ 

cLMz = ^L 2Re \L - Jo Ti (27F /(M2 - tl) - /o(6) K J" (21) 

Instead of the variable z used in the equations for p^-broadening it is equivalent to use the 

time variable t. In the large N c limit the coupling of the quark emitting a gluon is given by 

a a C F ^ a B N c 

TT 2 — 27T 2 ' 

The factor dt/X = podt counts the number of scatterers in the medium. The factor 1/L 
appears in ( pi]) because the spectrum is given per unit length. A is the mean free path of 
the quark and k = The k limits indicated in (|2T| ) eliminate the medium independent 
factorization contribution. 

These characteristic properties are now taken into account to allow the natural general- 
ization to the expanding medium. By properly specifying the time dependences we rewrite 
(HI) as 



udl a s N c 1 f r o+L dt 2 f** dti 
2Re ~ 1 



dudz 2ir 2 L 
d 2 b 



/(6;t 2j ti)-/ *(6 s t 2 ) 



A(t 2 ) J ro X(h) 

ui=oo(k=0) 



(22) 



where we assume that the quark hits the medium at time ro and travels a path length L, as 
discussed in the Introduction. The gluon propagation from ti — > t 2 is controlled by a Green 
function 



(23) 



The initial condition is 



/(&;ti,*i) = /„(Mi), 



(24) 



which is given by fl20|), where now fi = fi(ti). With the definition of the coefficient q(t) given 
in (|12D the emission spectrum is expressed in a rather symmetric form with respect to t\ 
and t 2 , namely by 

LodI 



Ujdl 1 { f T 0+L ft2 

^-j^-j = a s N c — i?e < y rft 2 y dh q{t 2 )q{t x ) 
r d 2 b r d 2 b' r r . r . , 

x / — J —b-b'G^^b^h) 



(25) 



In the logarithmic approximation the amplitude /(&; £ 2 , ti), and therefore the Green function, 
satisfies a Schrodinger equation for the 2-dimensional harmonic oscillator (actually with 
imaginary potential) MM. For fixed ti the equation reads 



f(b;t 2 ,tt), 



(26) 



with uJ^it) = iq{t)/uj. For an expanding medium the frequency of the oscillator u> (t) is 
time-dependent. In the Bjorken model |10| the temperature of the hot medium scales with 
time, as given in (0), which translates to 



"o(t) = ^oOo) 



TO 



(27) 



The explicit expression for the Green function is derived in Appendix B ||12|| , and given by 
eq.( |B13| ). In order to perform the 6-space integrations in (|25| ) it is convenient to change 
variables 



:,./■ = 2iisuJo(T )T ( — 
To 



l/2u 



(28) 



with the index v — 1/(2 — a), such that \ < v < 1. The 6-space integral is given by 
/ = 



d?b rd 2 b' r ,r r, > 

2 



2ut 



uj(2iiyu; (To)To) 2u 



2v-2 



l2 ' 



(29) 



in terms of modified Bessel functions I v (z) and K v (z) |13|. Inserting / and the time de- 
pendence of the coefficient q(t) as specified in (0) into the spectrum (|25|) a rather simple 
expression is obtained 



ujdl 
duodz 



x Re 



a s N c 1 1 
7T L 4z/ 2 

to+l dt 2 /-* 2 dti 

o h 



1 



*i [4_i(2;i)^_i(2; / ) - I v -i(zf)K v -i(zi 



(30) 



If we set 



XiJ = To 



^2 
T 



l/2f 



(31) 



and express the function K v (z) in terms of I± v (z), (excluding the case v = 1), we arrive at 

,i2 



cud I 

duodz 



a s N c 1 



7T 



2 sin7r(z/ — 1, 



7T 



f o g?x,- r 3 * dx 



1 



/to Xi Jt x f \I v _ l (2ii>uj x i )Ii- v (2ii>uj Xf) - (x { <-> a;/)]' 
where we put u; = c<j (r ) and f = r (l + L/r ) 1 ^ 2u for shorter notation. 

In order to compare with our previous result for the non-expanding plasma 
v = 1/2. The induced spectrum (^) then becomes 
ujdl a s N c „ f L+T " dxi r x * dxf 



(32) 



we take 



duodz 



Re 



Xi Jtq Xf 



Ii /2 {iuoX f )I^i /2 {iuj G x, i ) - (xi x/) 



(33) 
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Using h/ 2 (z) = ^-smhz/^/z and I -1/2(2) = ^-coshz/^/z |TJ gives 

1 2 



cudl a s N c r L r*i 

tie I dxj I dx 



duodz irL 



f 



ILOq 



sinh(iuJo(x j — Xi)) 



(34) 



where we have put r = 0. The remaining integrals can be performed explicitly, 



udl a.N r 



dudz 71L Jo Xj 



lU Q Xi 



a s N c 
71L 

OsNc 

irL 



Reiln 



tanh(zu; Xj) 
smh(iuj L) \ 1 
iu L I J 



(35) 



In 



sin(o;oL) 



u L 



This is the radiation spectrum in the iV c — > 00 limit, derived and discussed in 
hard quark entering the static QCD medium and radiating a soft gluon. 

On can easily go beyond the large N c limit and the soft gluon approximation in (p^ 
([35D. For a particle in an arbitrary colour representation R, should be replaced by 



^n^T^ll -X + —X 



for a 



and 



(36) 



where x is the gluon energy fraction x = u/E. In addition the r.h.s. of eqs.(|32|) and (|3~5l 
should be multiplied by 

2C R ,_ x 1 



1-X + 



(37) 



for a spin \ fermion, and by 



2C R (l + x 4 + (l-x) 



2(1 -x) 



(38) 



for a spin 1 particle (e.g. gluon, Cr = N c ) 



IV. ENERGY LOSS IN AN EXPANDING MEDIUM 



Next we integrate the radiation spectrum eq. (|32|) with respect to the gluon energy to in 
order to obtain the energy loss per unit length 
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dE _ r E ^ udl 
dz Jo dujdz ' 



(39) 



where we extend the limit E 
integration variables x and z 



oo. In analogy with the static case we introduce new 



2wuiQ{jQ)xi — i(l + i)x = x, z 



(40) 



leading to 



CO 



2v 2 q{r ) 



(41) 



Taking r = and performing the ^-integration 



J dxiXiq(T ) = ir 2 g(r ) ) = ~q{L)L 



(42) 



the energy loss can be written as 



dE _ 2a s N c 
dz 



71 



2z/sin7r(z/ — 1) 



7T 



q(L)L 



x Re 



1 dz r 00 dx 



x=0 



(43) 



/O Z Jo X 3 [Iv^fyh-vijbz) ~ I v -i(xz)h_ u (x) 

In order to obtain the subtraction term we expand the modified Bessel functions around 
x = [TB|], I v (z) — (hz) /r(^ + 1). This enables us to write (f4"5|) as 

2a s N c A/r , r r2r(z/+l)r(2-z/)sin7r(z/-l)l 2 
= — f_£ g (L )L — i — i - 

dz 



71 



7T 



dz 



o z \z v - 1 - z 1 -"] 2 



(44) 



where the function I(u, z) is defined in eq. (|Cl|) and evaluated in Appendix C. We integrate 
over the z variable (see Appendix C) and obtain the analytic expression for the energy loss 

"2r(i/ + l)r(2 - v) shi7r(z/ - l) n ~ 



dE a s N c 



7T 



(45) 



where the function 



IM 



4(1 -vf (2 -v) 



(46) 
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for | < v < 1 is derived in eqs.( |C4|) -( |C10|) . Notice that for v = 1/2 one recovers the energy 



loss for a quark traversing a static medium of size L, as discussed in 



dE 

dz 



static 



^mi*. (47) 



Eqs.(|4"5|) and ( fj7|) require that the tu-integration in (|39| ) be dominated by small x gluons. 
These formulas remain true beyond the large N c limit. The colour properties of the traversing 
particle are contained in the (transport) coefficient q(L) given in eq. (|l2"D . 
Using (flop and (|47|) one finds 

dE Qv 2 f dE \ 



dz 2 — v \ dz 



static , 



6 f dE 



(2 - a) (3 - 2a) \ dz 



a = 2--. (48) 



static / <s 



In case the quark is produced in the medium, rather than outside, 



— -2v (-— 

dz \ dz 



static 

dE 



2 — a \ dz 

replaces (ESI), where — ^ 



static / 



(49) 



also corresponds to a quark produced in the medium, and it 

static 

is 3 times the expression given in Q4"7D Q]. 

We notice that the limit u = 1 (a = 1) for an expanding ideal relativistic plasma can 
be taken. In this limit the maximal loss is achieved. It is bigger by a factor 6 for a quark 
produced outside ( 2 for inside ) than in the corresponding static case. In this comparison 
the temperature is taken after the expansion. The coefficient q(L) = q(T(L)) has to be 
evaluated at the temperature T(L) the quark finally "feels" after having passed the distance 
L through the medium which during this propagation cools down to T(L). It is remarkable 



that the initial temperature T of the hot medium does not enter the formulae ( pf ) and (f!9"D, 
although T is actually diverging in the limit r — > 0. 

As a consequence it is straightforward to generalize the relationship between energy 
loss and p^-broadening derived in |3[ for static nuclear matter to the case of an expanding 
plasma. We derive the relationships for a quark approaching the medium 
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dE a.N n 1 d 2 

L 7rrP±w( L )> (50) 



2 (2-ot)(3-2a) dL 

and for a quark produced in the medium 

a s N c 1 <9 2 
"^ = — 2^% P ^ (L) ' (51) 

relating the energy loss per unit distance in a hot expanding medium with the typical 
transverse momentum squared ( WM a jet receives in traversing a length L of a longitudinally 
expanding plasma. For a = the results of ||] are reproduced. 

APPENDIX A: PROPERTIES OF AN EXPANDING PLASMA 

Here we recall and briefly summarize the main properties of the space-time evolution 
of a hadronic fluid, which is produced by highly relativistic nucleus-nucleus collisions. We 
consider a hydrodynamical model and follow Bjorken [JlO ] in assuming one- dimensional lon- 
gitudinal expansion. 



In order to obtain the dependence of the fluid's temperature T = T(r) on the proper-time 
r we use the conservation law for the entropy density s, 

ds/dr + s/t = 0. (Al) 
We take into account the thermodynamic equation for the pressure 

dp/dT=s(T(r)), (A2) 
and express p in terms of a monotonically increasing function of temperature, n(T), 

P = (T)T 4 . (A3) 



Defining the parameter 



T dn(T) fk . 

Ai 3 MTj^r- (A4) 
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which we assume to be temperature independent^. It follows from eqs. ([Al| -| 

dJ ^ = -^ vl = {Z + ^l\ (AS) 



with v s the sound velocity. In the approximation of v s = const eq. (|A5|) gives 



(T/T ) 3 = (r /r) a , witha = 3v 2 . (A6) 

The parameter a is bounded by < a < 1, where a = means constant temperature and 
a static medium, a = 1 is an ideal relativistic plasma. 

In perturbative thermal QCD [|ll|] the parameter Ai turns out to be small, indicating 
small deviations from ideal gas behaviour. For the case of a gluon gas 



165 /a s \ 2 



(^) (l + 0(v^)), (A7) 



in terms of the QCD coupling constant a s , which at very high temperatures should be 
evaluated at the scale T. 

APPENDIX B: GREEN FUNCTION 



In order to discuss the solution of eq.(|26|) we make a logarithmic approximation and 
assume that ln(l/6 2 ) is slowly varying for small b 2 . We then have to solve the Schrodinger 
equation for a two-dimensional harmonic oscillator with time-dependent frequency. Using 
variables familiar from quantum mechanics the equation is 



f(b,t), (Bl) 



where the mass of the oscillator is identified with the energy of the emitted gluon, m = — oo, 
and the time dependence of the frequency is given by the power behaviour 

u 2 (t)=u 2 (t )(t /ty, (B2) 



2 Possible T dependences of Ai are sketched in ref. [10 
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where the parameter a is discussed in Appendix A, and 

iq(t ) 



^o(*o) 



(B3) 



The Green function of the Schrodinger equation ( |B1|) can be written in the form |L2 



G(b,t; b',t' 



m 



— exp {iS c i(b,t;b',t')} , 



2iriD(t,t' 

where the fluctuation determinant satisfies the homogeneous differential equation 

d 2 



(B4) 



dt 2 



D{t,t')+u 2 {t)D{t,t') = 0, 



(B5) 



with the initial conditions 



D(t',t')=0, ^-D(t,t')\ t=t , = 1. 
dt 



(B6) 



The classical action 5* c i in (|B4| ) is determined by the classical path b c \(t) obeying 

d 2 



dt 2 



b cl (t)+LU 2 (t)b cl (t)=0, 



(B7) 



with 



b cl (t) = b and b d (t') = b' 



(B8) 



It follows that 



m 



S cl (b,t;b',t') = - 



bd(t) ■ jb d {t) 



(B9) 



The explicit solution of (|B5| ) with (p^) is found in terms of modified Bessel functions I v {z) 
and Kjz) ill to be 



D(t,t') 



2ut 



[2tvuj (t )t } 2 » 
where we introduce the variables 



[zz'Y \I v {z)K v {z') - K v {z)I v {z')\ 



(BIO) 



z = z(t) = 2iuou {t )t {t/t ) 1/2 \ z' = z(t') 



(Bll) 



with 
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i/= 1/(2 -a), 



(B12) 



such that 1/2 < v < 1. 

Using the solution to (|B7|) and ( |B~8D in ( [B9D gives the Green function (|B4|) as 



— 10J 



with the coefficients 

Cl = z^/z)" [I v ^{z)K v {z') + /^(z)/^')] , 
c 2 = z'(z/z'Y [K v {z)I v ^(z') + /^^(z')] . 



(B13) 



(B14) 



The case v = 1/2 is especially easy to handle, and allows a direct comparison with the 
results already obtained in H||. The variables given in ( |B11|) become z(z') = iuot(t') with 
lu = Jifi 2 v/Xuj. The functions h/2(z) and i^i/ 2 (^) are expressed in terms of hyperbolic 



functions |T3||, so that the determinant ( |B10| ) simplifies to 



D(t, t') = — sinh (z — z ) = — sin uo(t — t') 



We note that the Green function ( |B13| ) is time-translational invariant for u — 1/2. It 
correctly reproduces the result of eq.(5.6) in @. 



APPENDIX C: THE INTEGRALS z) AND I{y) 



Here we evaluate the integral 

f°° dx { 
I{u,*)=Ref ^ 1 



[r»r(2 - „)]- 



n2 [ ' 



(Cl) 



[h-„{x)I v -i{xz) - h- u {xz)I u -\{x)Y 
with x = z(l + by using the integration contour Ci r ..4 in the complex z = (x + iy)- 
plane, which we already introduced in |2| (see Appendix D), and which for convenience is 
reproduced here in Fig. 1. 
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FIGURES 



y = Imz 




x = Rez 



FIG. 1. Integration contour for the integral l{y,z) 



Following the detailed discussion given in |2| the integral is performed by calculating the 
residue of the pole at x = 0. Since the contribution along C 2 vanishes, we find the result for 
I(u, z) by adding the contributions from the paths C\ and C3, i.e. 



H Residue 



.1-' 



■{•/•} 



|z=0; 



(C2) 



leading to 



7T2; 



si;) (^- 2 - * 2 -") + J - *~ 



(C3) 



For the reader who may not be convinced by the above arguments we note that we have 
evaluated the integral ( |C1|) numerically using the program Mathematica |14| . Stable results 
were obtained agreeing with ( |C3|) for a large domain of z and u, 0.1 < z < 0.8 and | < v < 
0.95. 

In eq.(pE5|) we stated that the energy loss —dE/dz is proportional to the integral I(y) 
defined by 



IT Jo Z \ Z V - X - Z l ~ u ] 



dz 



(2-iz) 



\Z V - 2 - Z 2 ' U \ + i \ Z " - Z~ V \ 



yU—1 ~1 — 1/1 



(C4) 
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For the special case v = 1/2 



7(1/2) = ~ C dz = 2/3. 
6 Jo 



(C5) 



For z/ in the interval \ < v < 1 we evaluate 7(z/) as follows. We change the integration 
variable to 



t = z 



2(1-") 



(C6) 



We regularize the integrand near t = 1 by (1 — t) 3 — ► (1 — t) 3+E , £ > 0, and arrive at 



/e(z/) 2v{l-v){2-v) Jo 



i z/ 



+ {2-u) 



t~ - t 



(C7) 



(l-t) 3 - £ 

where the limit e — > is to be taken after the integration. Using the Euler beta-function 
gives 

1 { v 2- is 



2v{\ -v){2-v) \e-2 (e-l)(e-2) 

z/(2 - u)T(e + 1) 
(l-i/)Me-l)(e-2) 



(C8) 



One can easily check that I e {y) is regular at e — 0. Using [TT3 



r(z)/r(e + z) — > 1 - e*l>{z) + 0(e 2 ) 



in terms of the digamma function ip(z), and with the recurrence formula ip(z + 1) 
1/z, we finally obtain 

Uu) — >Uv) = — 1- -, -<V<1. 

£K e^o y 1 4(1 -z/) 2 (2-z/)' 2 - 



(C9) 
(CIO) 



This is in agreement with ( |C5|) for v — 1/2 
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